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Abstract We present a set of Bell inequalities for multi-qubit quantum systems. These Bell inequal- 
ities are shown to be able to detect multi-qubit entanglement better than previous Bell inequalities 
such as WWZB ones. Computable formulas are presented for calculating the maximal violations of 
O | these Bell inequalities for any multi-qubit states. 
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Based on the Einstein, Podolsky and Rosen (EPR) gedanken experiment, Bell presented an in- 
. equality obeyed by any local hidden variable theory [lj. It turns out that this inequality and its 

generalized forms are satisfied by all separable quantum states, but may be violated by pure en- 
tangled states and some mixed quantum states (2j 12111 13 0|71. Thus Bell inequalities are of great 
importance both in understanding the conceptual foundations of quantum theory and in investigating 
t— I " quantum entanglement. Bell inequalities are also closely related to certain tasks in quantum infor- 

^ | mation processing, such as building quantum protocols to decrease communication complexity [8] 

and providing secure quantum communication [9]. Due to their significance, Bell inequalities have 
been generalized from two qubit case, such as the Clauser-Horne-Shimony-Holt (CHSH) inequal- 
ity ifTOl to the N-qubit case, such as the Mermin-Ardehali-Behnskii-Klyshko (MABK) inequality 
ifTTl [T2Tl . and to arbitrary d-dimensional (qudit) systems such as the Collins-Gisin-Linden-Masser- 
Popescu (CGLMP) inequality |fl3l . However, except for some special cases such as bipartite pure 
states (2j[3j|6l, three-qubit pure states Bl [T4l . and general two-qubit quantum states [15 ], there are 
no Bell inequalities yet that can be violated by all the entangled quantum states, although it is shown 
recently that any entangled multipartite pure states should violate a Bell inequality Q. Thus it is of 
great importance to find more effective Bell type inequalities to detect the quantum entanglement. 

In this paper, we study Bell inequalities for both pure and mixed multi-qubit systems. We 
propose a series of Bell inequalities for any Af-qubit states (N > 3), and derive the formulas of the 
maximal violations of these Bell inequalities. This gives a sufficient and necessary condition which 
is also practical for any multipartite qubits quantum states. It is shown that the Bell inequalities 
constructed in this paper are independent of the WWZB inequality and Chen's Bell inequalities 
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constructed in (0, i.e. they can detect some entangled states which fulfill both the WWZB inequality 
and Chen's Bell inequalities. 

Consider an A-qubit quantum system and allow each part to choose independently between 
two dichotomic observables A,-, A. for the zth observer, i = 1,2, ...,N. Each measurement has 
two possible outcomes 1 and -1. Quantum mechanically these observables can be expressed as 
A; = dj(f, A'- = a t cr, where cr = (cry, cr 2 , C3) are the Pauli matrices and a,, a i are unit vectors, 
i = 1,2,- •• ,N. 

The CHSH Bell inequality for two-qubit systems is given by 

\(B 2 )\ < 1, (1) 

where the Bell operator B 2 = 5[AiA 2 + A\A 2 + A\A' 2 - A'jAj]. In HH Horodeckis have derived 
an elegant formula which serves as a necessary and sufficient condition for violating the CHSH 
inequality by an arbitrary mixed two qubits state. 
The WWZB Bell operator is defined by 

B n WZB = ^n Z S( Sl ,s 2 ,---,s N ) ^ ^•••4"<=i°#A ( 2 ) 

si,S2,--- ,sn=±\ k\,ki,--- ,kf/=±l 

where S(s\,s 2 , - ■ ■ , sn) is an arbitrary function of s;(= ±1), i = 1, ...,N, taking values ±1, 0/(1) - 
Aj and 0/(2) = A\ with kj = 1, 2. It is shown in HTHHl that local realism requires \(B N )\ < 1. The 

MABK inequality is recovered by taking S(si,s 2 ,--- , s^) = V2cos[(^i + s 2 H 1- - N+ 1)/|] 

in (f2]). In ifTTl |T8l the authors also derived a necessary and sufficient condition of violation of this 
inequality for an arbitrary N-qubit mixed state, generalizing two-qubit results in lfl6l . However, 
when using the results to obtain the maximal violation of the WWZB inequality, one has to select 
a proper set of local coordinate systems and a proper set of unit vectors, which makes the approach 
less operational. 

Employing an inductive method from the (N- l)-partite WWZB Bell inequality to the A^-partite 
inequality, a family of Bell inequalities was presented in Ifl2l . The Bell operator is defined by 

B N - B™ ® j(A N + A N ) + ® 1 -{A N - A N ), (3) 

where B^ W ^ ZB represents the normal WWZB Bell operators defined in (f2]), Zjv-i be the identity 
operators acting on first (A^- 1) qubits. Such Bell operators yield the violation of the Bell inequality 
for the generalized GHZ states, = cos a|00 • • • 0) + sin a\ 1 1 • • • 1 ), in the whole parameter region 
of a such that cos a + and sin a + 0, and for any number of qubits, thus overcoming the drawback 
of the WWZB inequality. In the three-qubit case, one can construct three different Bell operators 
from B 2 by using the approach of Q. The corresponding three Bell inequalities can distinguish 
full separability, detailed partial separability and true entanglement |fl9l . However, the maximal 
violation of this Bell inequality is unknown for a generally given three-qubit state. 

We start with constructing a set of new Bell inequalities for any A^-qubit quantum systems by 
iteration. First consider the case Af = 3. As a two-qubit CHSH Bell operator S 2 can act on two of 
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the three qubits in three different ways, we can have three Bell operators, 

si = (s 2 y ® \(Ai + a;.) + {hi ® l(Ai - a;.), / - 1, 2, 3, (4) 

where (B2) 1 and (h)' are the two-qubit CHSH Bell operator and the identity operator acting on the 
two qubits except for the ith one. For N > 4, the Bell operators can be similarly obtained, 

with i = 1, 2, • • • ,iV and 7 = 1, 2, • • • , ^~^ ! . Here (S'Lj)' denotes the jth Bell operator acting on 
the (N - 1) qubits except for the ith one. (In-i)' stands for the identity operator acting on the (N - 1) 
qubits except for the jth one. There are totally ^ Bell operators. 

Theorem 1: If a local realistic description is assumed, the following inequalities must hold, 

\(S k N )\ < 1, (6) 

where k e {1,2,- •• , f }. 

Proof: We prove the theorem by induction. Note that for two qubits systems, local realism 
requires that K#2)l < 1 as shown in £[]). Assume that a local realistic model has lead to \{S k N l )\ < 1 
with k € {1, 2, • • • , -^p- }■ We consider the N-partite systems. If A; and A . are specified by some 
local parameters each having two possible outcomes -1 and 1, one has either |A,- + A,| = 2 and 
|A; - A.| = 0, or vice versa. For any k e {1,2, • • • , ^ j, from ((5]) we have that 

l<s*>l - \(st l) 2 j >i = K(^_o ! '®-(A£+a ! .) + (W®-(a ! -a'.))| 



< k^.j)')! ® \\aAi + a;.)>i + i|<cA, - a;.)>i < 1. 

□ 

It is shown in EUll that violation of a Bell inequality gives rise formally to a kind of entanglement 
witness. Moreover, the separability criterion and the existence of a description of the state by 
a local hidden variable theory will become equivalent when one restricts the set of local hidden 
variable theories to the domain of quantum mechanics. Thus one can use the Bell inequalities as the 
separability criteria to detect quantum entanglement. We remark that any A^-qubit fully separable 
states also satisfy the inequality ©. For N > 4, the operator & N _^ is derived from B' N _ 2 - Thus 
S' N are different from the Bell operators in [12] where S' N _ l is the Bell operator in the WWZB 
inequality. The following example will show that our Bell inequalities in (O are independent from 
the WWZB inequalities and that in Ifl2l . and our new Bell inequalities can detect entanglement 
better than they can. 

Example Consider a four-qubit pure state \if/) = \<p) ® |0>, where \<f>) - cosa|000) + sin a\ 111), 
a e [0, ■jg]. It has been proved ED that for sin 2a < \ (i.e. a e [0, ^]), the WWZB Bell inequalities 
cannot be violated by the generalized GHZ state \<p). According to the result in EH . the WWZB 
inequalities operator B^ WZB and the Bell operator B4 in Ifl2l satisfy the following relations, 

\{wr ZB \n * MBV^m * h (7) 
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MBaW>\ * MB« W£B \4»\ < 1. (8) 

Therefore both the WWZB Bell inequalities and the inequalities in [12] can not detect entanglement 
of ItA). 

Nevertheless the mean values of the Bell operator & 2 in © is V2 sin 2 2a + cos 2 2a which is 
always larger than 1 as long as \<p) is not separable. Therefore the entanglement is detected by our 
Bell inequality ©. 

To identify the nonTocal properties of a quantum state, it is important to find the necessary and 
sufficient conditions for the quantum state to violate the Bell inequality. Now we investigate the 
maximal violation of the Bell inequalities Q. We first consider the N = 3 case. In this situation, 
© gives three operators, 

S\ = (S 2 ) 1 ®i(A 1 +A' 1 ) + (/ 2 ) 1 ®i(A 1 -A' 1 ), (9) 
®l = (&2) 2 ®j(A2+A 2 ) + (I 2 ) 2 ®^(A 2 -A 2 ), (10) 

£3 = (&2) 3 ®\(A3+A 3 ) + (I 2 ) 3 ®^(A 3 -A 3 ), (11) 

where {S 2 f = \(A 2 A 3 + A' 2 A 3 + A 2 A 3 - A' 2 A' 3 ), (S 2 ) 2 = \{A { A 3 + A[A 3 + AjA 3 - A\a' 3 ) and 
(S 2 ) 3 = ^(AiA 2 + A\A 2 + A1A2 - AjAl). Let p be a general three-qubit state, 

3 

P = 2j TijUO-iCTjCTk, (12) 
i,#=0 

where cr = 7 2 is the 2x2 identity matrix, o~ ( - are the Pauli matrices, and 

Tijk = ^Tripo-jO-jO-k). (13) 

Theorem 2: The maximum of the mean values of the Bell operators in (O, (flOl and (fTTT ) satisfy 
the following relations, 

max 1(^)1 = 8maxM}(^) + 4(^) + ll^ || 2 -<^,f 1 > 2 }i (14) 
max|(S 2 >| = 8maxU 2 (^ 3 ) + + lir 2 || 2 - (4 f 2 > 2 }?, (15) 

max|(^>| = 8max{^(£ 3 ) + A \(b 3 ) + ||7* || 2 - (b 3 , f 3 > 2 }5, (16) 

where (.,.) denotes the inner product of two vectors, ||x|| stands for the norm of vector x. The 
maximums on the right of (fT4l ). ([TBI and ([TBI are taken over all the unit vectors b 3 . Given a three- 
qubit state p, one can compute by using the formula in ( [TBI . Then ^(63) and A' 2 (b 3 ) are defined 
to be the two greater eigenvalues of the matrix MjMi with Mj - YX=\ b 3 T' k , i = 1, 2, 3, with respect 
to the three Bell operators in ©, (flOl and (fTTb . Here Tl, I = 1,2, 3, are matrices with entries given 
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b y ( T k^ij ~ Tki P ( T k^j - T &j and ( T pij = T Uk- T™ y m = 1 ' 2 ' 3 we defined to be vectors with 
entries (7^,)* - T m , (f^) k - T oko and (7^)* = Too*. 

Proof: We take (fTTT) as an example to show how to calculate the maximal violation. The max- 
imal violation for the Bell operators (O and (fTOl ) can be computed similarly. A direct computation 
shows that 

<B\ = i[(Ai + A\)A 2 + (A l - A\)A 2 ](A 3 + A 3 ) + {l 2 f ® I(A 3 - A 3 ) 

= i(Ai + A\)A 2 (A 3 + A 3 ) + i(Ai - A\)A 2 (A 3 + A' 3 ) + (7 2 ) 3 ® ^(A 3 - A^ 

1 3 , , 1 3 

= 4 E ^1 + ( a 'i^ a 2^4 + (a' 3 ) k ]o-iO-jO- k + - 2_j K - (a'i)'](a2) 7 [ fl 3 + i^t^i^j^k 

i,j,k=l i,j,k=l 
1 3 

+ -^[4-(fl 3 )']/4® ( 7 t . (17) 
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For any given unit vectors a\ and a\, there always exist a pair of unit and mutually orthogonal 
vectors fti, ft 1 and e [0, |] such that 

-* -> -> -> 

ci\ + a'l = 2cos 0b\, <?i — a 1 = 2sin0ft'i. (18) 

^» 

Similarly for a* 3 and a'3, we have 

a 3 + a 3 = 2 cos (f>b 3 , a 3 - a 3 = 2 sin <pb 3, (19) 

where ^3 and ft' 3 are orthogonal vectors with unit norm and cp e [0, |]. 

By inserting (fT2l) into (1221) we get the mean value of the Bell operator (fTTT) . 



<^> = Tr(pB\) 
1 3 
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+ I E [fl i - («i)'K«2) 7 '[«3 + (^AWK^**) 

1 3 

+ 7^[«3" {a' 3 ) k ]T m Tr(lj ® o\) 



k=\ 

3 3 



= 8 b\b 3 a J 2 T ijk cos 6>cos + 8 ^ {b x ib\(a 2 yTi jk sin 0cos + 4 ^(ft^Too/t sin <f>. 

i,j,k=l i,j,k=l k=l 

Let T^,k = 1,2,3, be the matrux with entries given by (T^)ij = Tij k and T^ Q a vector with 
components (f^k = ^oo/t- The maximal mean value of the Bell operator (fTTT> can be written as 

3 3 
max<S^) = 8max[<ft^,^ft^r|<J 2 >cos0cos0 + {b v ^b k 3 Tld 2 ) sin 9 cos (p + (b 3 ,f^ Q ) sin (f>] 

k=l k=l 
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8max{[{bi,^b k 3 T 3 d 2 )cos9 + {^j,^ ^7^4) sin 0] 2 + (b' 3 ,T 3 ) 2 }? 
k=i k=i 

3 3 

8max{[<£j, j^r t 3 a 2 >cos0 + (b\,^b k 3 T 3 k a' 2 ) sinO] 2 + ||f 3 || 2 - (h,?^) 
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&=1 fc=l 

= 8 max{i 3 (£ 3 ) + A 3 2 (b 3 ) + ||f 3 || 2 - (b 3 , f 3 ^) 2 }?, (20) 

which proves (fT6l ). In (1201) we have used the fact that the maximum of x cos + y sin taking over 
all is a/jc 2 + y 2 . Formulae (fl4l) and (fT5T > can be similarly proven. □ 

Remark: According to the symmetry of the operator S 3 , the equation (|20l also provides the 
minimum of the operator (fTTT >. achieved by -S 3 . 

Since b 3 is a three dimensional real unit vector, one can always calculate the exact value of the 
maximum for any given three qubits quantum state. For example, for the generalized three-qubit 
GHZ state, \GHZ) - cos a|000) + sin a\ 1 1 1), by selecting some proper direction of the measurement 
operators, i.e. a,s and (a ),s, the maximal mean value of the Bell operator in (TTTb is shown to be 
sin 2a + cos 2 2a. From our formulae in Theorem 2 one can show that the result are in 
accord with that in [12]. For three-qubit W state, \W) = -^=(|100> + |010> + |001», our mean value 
is 1.202, which is also in agreement with that in [12J. However, our method can be also used to 
calculate the mean value of the Bell operators in ©, (flOl and (TTTb for any three qubits quantum 
states. For instance, we consider the mixture of \W) and \GHZ), 

p = X \W)(W\ + (1 - x)\GHZ)(GHZ\, (21) 

where < x < 1. We have the maximal mean value of the Bell operator in (fTTIl.see FigQ] where 
f(x) stands for the maximal mean value of the operator (fTTT > for the mixed state p. For < x < 0.33 
and 0.82 < x < 1, f(x) > 1 and p is detected to be entangled. 

For four-qubit systems, the Bell operators ([5]) have four different forms. We take S 12 as an 
example to investigate the maximal violation of the corresponding Bell inequality. Note that 

Sf = (& 3 ) 4 ®^(A 4 +A 4 ) + (I 3 ) 4 ®^(A 4 -A 4 ) 

1 , 1 . , 

- -(Ai + A l )A 2 (A 3 +A3XA4 + A 4 ) + -(Aj - Aj)A 2 (A 3 +A3XA4 +A 4 ) 

o 

1 , 1 

+-h ® (A 3 - A 3 )(A 4 + A 4 ) + -I 6 <g> (A 4 - A 4 ) 

1 3 

= 7 ^ [a\+(a\)']a J 2 [a 3 + (a' 3 ) k ][a[ + (a'J](T i ajcr k o-i 
Uj,k,i=i 

1 3 

1 3 , , 1 3 



fc,Z=l /=1 
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f (X) 




Figure 1: The maximal mean value of the operator (fTTT ) for the mixed state p in (|2TT >. f(x) stands 
for the maximal mean value and x is the parameter in p. 

Let p be a general four-qubit quantum state, 

3 

P = 2j T ijklO~iO~jO~kO-l, (23) 

iJ,*,/=o 

with r ; yjy = r{pcr iCTjcrkO-i)- The mean value of S l 4 2 can be derived by the following deduction. 
1 3 

(S\ 2 ) = - J] K + («i) ! ']«2[«3 + (4)*]t°4 + (aJmjuTiiajajaiaj) 

+ 8 2 [a i ~ ("l)'"]^^ + (« 3 ) fc ][«4 + (a^mjuTriajajaiaj) 

ij,kj=l 

+ 4 Z [a 3 ~ (4)*] [4 + («4)']rooHrr(/4 ® c^of) 

j 3 

z=i 

3 

- 2 4 ^ ft j a^b^b^TijM cos cos a 3 cos a 4 

3 

cos »3 cos a\ 

3 3 

+2 4 2] (b' 3 ) k b' 4 T m i sin or 3 cos a 4 + 2 4 ^(ft 4 )'r oo/ sin a 4 , 
k,l=i l=l 

where we have assumed that a,- + a', = 2 cos a,- ft;, a,- - a',- = 2 sin ar, ft',, or, € [0, |]. 
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The maximum of the mean value can be derived to be 

3 

max(!B\ 2 ) = 2 4 max[(b\ , ti^b^T^fy) cos a \ cos #3 cos ®4 

k,i=i 

3 

+{b x , y b 3 b 4 T kl a 2 ) sin ct \ cos (23 cos a 4 

+(& 3 , ^00^4) sin of 3 cos ^4] + (fc 4 , Tqqq) sin ar 4 ] 
3 3 
= 2 4 max{<^, J] b\b l jlfa 2 ) 2 + (b\, £ b\b l A T k l2 ct 2 ) 2 + <£,, r o %> 2 + <£ 4 , f,J 2 > 2 }5 

- 2 4 maxM} 2 (£ 3 £ 4 ) + 4W 4 ) + ll^ 2 ^!! 2 - (h, T X2 b 4 ) 2 + \\f£ \\ 2 - $ 4 , f^) 2 }K 

where A l2 (bi,b 4 ) and A l2 (bT,b 4 ) are the two greater eigenvalues of the matrix (M 12 )^M 12 , M n = 
Y<1 i = i b^b^T^; T k f stand for the matrices with enuies (Tfflij = Tiju with i, j, k,l - 1, 2, 3; T^j is a 
matrix with entries (T^) k i = 7oo/y, and 7^ is a vector with components (T^ Q )i = 7ooo/, I - 1, 2, 3. 
The maximum in the last equation is taken over all the unit vectors Z? 3 and b 4 . 

In terms of the analysis above, for four-qubit systems we have the following Theorem. 

Theorem 3: The maximum of the mean values of the Bell operators in (f5]) for four qubits 
systems are given by the following formula: 

max |<fl*>| = 2 4 maxU'f (£ 3 £ 4 ) + %$Ja) + \\T^ 4 \\ 2 - T™b*4? + ||7^ || 2 - ($ 4 , f^fp. (24) 

The maximums on the right side are taken over all the unit vectors Z? 3 and b 4 . Here A™(b 3 b 4 ) 
and A 2 {t 3 $ 4 ) are the two greater eigenvalues of the matrix {M m )^M m , M m = £j^ = i b^b^T™, m = 
1, 2, • • • , ^; are the matrices with entries (T^j = T Uj , (r 2 )y = T ukj , (7*)y = T lijk , (r 4 ) i; - = 
Tkiij, (Tfyj = T nkj , {Tfyij = T iljk , {T 7 kl )ij = T mj , (7*) y = T Wj , {T^ij = T iflk , (7^% - 7*yz, 
(Ty)ij = Tj k ji and (Tj£)ij - Tij k u i,j = 1,2,3 and / = 0,1,2,3; 7^ stand for the vectors 

with components (T' m ) x = T^oo, ( T ooo)x - T 0x00 , (T&x))* ~ To0x0 ' ( T oo<? x = T °° 0x ' 1 = 1?2 ' 3 ' 
j = 4,5,6, k = 7,8,9,/= 10, 11, 12 and % = 1,2,3. 

As an example, consider the 4-qubit W state \W) = ^(|1000> + |0100> + |0010> + |0001», by 
using the formula (1241) one gets the maximal mean value max \(S 12 )\ = 1.118. For the mixed state 
p = jgl + (1 - jc)|W)<W|, entanglement can be detected by (f24]) for < x < 0.106. 

Generally, for any N-qubit quantum state, the maximal mean values of the Bell operators in © 

can be calculated similarly by using our approach above. For example, the maximal mean value of 

m 

can be expressed as 

A" 

max \(Sj)\ = 2 N m^{Af(b 3 ■ ■ ■ b N ) + A%{$ 3 ■ ■ ■ b N ) + ||7V,vl| 2 - <^3, ?45-n) 2 

+ 1175-Arll 2 - <^4, 75-isr> 2 + • - • + IITVll 2 - <^v, 7V> 2 }^, (25) 

where A\(b 3 ■ ■ -b^) and ^2(^3 • • -fcjv) are the two greater eigenvalues of the matrix M^M, with 
(M)ij = ... ( - =1 £3 • • • b'^Tiji 3 „.j N the entries of matrix M; T45...N, Ts-n and 7% are vectors with 
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components (f 4 5-tf)* = 2L-,w=i b 4 ' ' ' b N T ^h,-,iN'(^5-N)k = , ^ 5 • • ' ^^000^,-,^ 

and (i at)j = Tooo-ofc respectively. The maximum on the right side is taken over all the unit vectors 
~* ~* ~* i — i 

03,04, •• • , Ojv- The other mean values of the Bell operators in ((5J) for Af-qubit states can be obtained 

-> 

similarly. By expressing the unit vectors as (cos 9\ cos fa, cos Ok sin fa, sin 9k), k - 3, • • • , iV, our 
formulas can be used to compute the maximal violation by searching for the maximum over all 9k 
and fa, either analytically or numerically. 

In conclusion, we have presented a series of Bell inequalities for multipartite qubits systems. 
These Bell inequalities are more effective in detecting the non-local properties of quantum states 
that can not be described by local realistic models. Formulas are derived to calculate the maximal 
violation of the Bell inequalities for any given multiqubit states, which gives rise to the sufficient 
and necessary condition for the violation of these Bell inequalities. For a fixed multiqubit state, 
one can optimize the mean values of the Bell operators over all measurement directions. Our Bell 
operators only involve two measurement settings per site, which meets the simplicity requirements 
of current linear optical experiments for nonlocality tests. Moreover, our formulas for the maximal 
violation of the Bell inequalities fit for both pure and mixed states, and can be used to improve the 
detection of multiqubit entanglement. 
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